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Dedicated underground experiments searcliing for dark matter have little sensitivity to GeV and 
sub-GeV masses of dark matter particles. We show that the decay of B mesons to K(K*) and 
missing energy in the final state can be an efficient probe of dark matter models in this mass range. 
We analyze the minimal scalar dark matter model to show that the width of the decay mode with 
two dark matter scalars B — > KSS may exceed the decay width in the Standard Model channel, 
B — > Kvv, by up to two orders of magnitude. Existing data from B physics experiments almost 
entirely exclude dark matter scalars with masses less than 1 GeV. Expected data from B factories 
will probe the range of dark matter masses up to 2 GeV. 



INTRODUCTION 



Although the existence of dark matter is firmly estab- 
lished through its gravitational interaction, the identity 
of dark matter remains a big mystery. Of special interest 
for particle physics are models of weakly interacting mas- 
sive particles (WIMPs), which have a number of attrac- 
tive features: well-understood mechanisms of ensuring 
the correct abundance through the annihilation at the 
freeze-out, milli-weak to weak strength of couplings to 
the "visible" sector of the Standard Model (SM), and as 
a consequence, distinct possibilities for WIMP detection. 
The main parameter governing the abundance today is 
WIMP annihilation cross section directly related to the 
dark matter abundance. In order to keep WIMP abun- 
dance equal or smaller than the observed dark matter 
energy density, the annihilation cross section has to sat- 
isfy the lower bound, (TannWrcl > 1 pb, (see e.g. 0). In 
all WIMP models studied to date the annihilation cross 
section is suppressed in the limit of very large or very 
small mass of a WIMP particle S. This confines the 
mass of a stable WIMP within a certain mass range, 
Wmin ^ < TTijjjaxj which wc refer to as the Lee- 
Weinberg window This window is model-dependent 
and typcally extends from a few GeV to a few TeV. If the 
neutralino is the lightest stable supersymmetric particle, 
mmin — 5 GeV H but in other models of dark matter 
Wmin can bc lowered 



Recently, WIMPs with masses in the GeV and sub- 
GeV range have been proposed as a solution to certain 
problems in astrophysics and cosmology. For example, 
sub-GeV WIMPs can produce a high yield of positrons 
in the products of WIMP annihilation near the centres of 
galaxies 0, which may account for 511 KeV photons ob- 
served recently in the emission from the Galactic bulge 
0- GeV-scale WIMPs are also preferred in models of 
self-interacting dark matter Q that can rectify the prob- 
lem with over-dense galactic centers predicted in numer- 
ical simulations with non-interacting cold dark matter. 



Dedicated underground experiments have little sensi- 
tivity to dark matter in the GeV and sub-GeV range. 
Direct detection of the nuclear recoil from the scatter- 
ing of such relatively light particles is very difficult be- 
cause of the rather low energy transfer to nuclei, Ai? ~ 
w^m|/mNuci Ss 0.1 KeV, which significantly weakens ex- 
perimental bounds on scattering rates below ms of few 
GeV, especially for heavy nuclei. Indirect detection via 
energetic neutrinos from the annihilation in the centre 
of the Sun/Earth is simply not possible in this mass 
range because of the absence of directionality. There- 
fore, the direct production of dark matter particles in 
particle physics experiments stands out as the most reli- 
able way of detecting WIMPs in the GeV and sub-GeV 
mass range. 

The purpose of this work is to prove that B decays can 
be an effective probe of dark matter near the lower edge 
of the Lee- Weinberg window. K decays can also be used 
for this purpose, but B decays have far greater reach, 
up to ms ^ 2.6 GeV. In particular, we show that pair 
production of WIMPs in the decays B ^ K{K*)SS can 
compete with the Standard Model mode B — > K{K*)vv. 
In what follows, we analyze in detail the "missing enere 
processes in the model of the singlet scalar WIMPs 
hd and use the existing data from B physics experiments 
to put important limits on the allowed mass range of 
scalar WIMPs. 

The main advantage of the singlet scalar model of dark 
matter is its simplicity, 



(1) 



where H is the SM Higgs field doublet, vew — 246 
GeV is the Higgs vacuum expectation value (vev) and 
h is the field corresponding to the physical Higgs, H — 
(0, {vew + hl)/V^)- It is important to recognize that the 
physical mass of the scalar S receives contributions from 
two terms, m|. = TOq + Xv^w^ ^^'^ '^^^ be small, even if 
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FIG. 1: Feynman diagrams which contribute to B meson 
decays with missing energy. 



each term is on the order ±0(w|;^). Although admit- 
tedly fine-tuned, the possibility of low ms is not a priori 
excluded and deserves further studies as it also leads to 
Higgs decays saturated by the invisible channel, h — > SS 
and suppression of all observable modes of Higgs decay at 
hadronic colliders 4] . The minimal scalar model is not a 
unique possibility for light dark matter, which can be in- 
troduced more naturally in other models. If for example, 
the dark matter scalar S couples to the Hd Higgs doublet 
in the two-Higgs modification of iQ), XS^Hj^Hd, the fine- 
tuning can be relaxed if the ratio of the two electroweak 
vevs, tan/3 = {Hy)/{Hd) is a large parameter. A well- 
motivated case of tan/3 ~ 50 corresponds to {Hd) ^ 5 
GeV, and only a modest degree of cancellation between 
TOq and X{Hd)'^ would be required to bring ms in the 
GeV range. More model-building possibilities open up 
if new particles, other than electroweak gauge bosons or 
Higgses, mediate the interaction between WIMPs and 
SM particles. If the mass scale of these new particles is 
smaller than the electroweak scale sub-GeV WIMPs 
are possible without fine-tuning. 



PAIR-PRODUCTION OF WIMPS IN B DECAYS 

The Higgs mass mh is heavy compared to ms of in- 
terest, which means that in all processes such as an- 
nihilation, pair production, and elastic scattering of S 
particles, A and mh will enter in the same combination, 
A^m^^. In what follows, we calculate the pair-production 
of S particles in B decays in terms of two parameters, 
\^ /mf^ and ms, and relate them using the dark matter 
abundance calculation, thus obtaining the definitive pre- 
diction for the signal as a function of ms alone. 

At the quark level the decays of the B meson with miss- 
ing energy correspond to the processes shown in Figure 
n The SM neutrino decay channel is shown in Figure 
^ and^. The 6 — s- s Higgs penguin transition. Figure 
produces a pair of scalar WIMPs S in the final state, 
which likewise leave a missing energy signal. In this sec- 
tion, we show that this additional amplitude generates 
6 — > sSS decays that can successfully compete with the 
SM neutrino channel. 

A loop-generated h — s— Higgs vertex at low momen- 
tum transfer can be easily calculated by differentiating 
the two-point b s amplitude over vew- We find that 



to leading order the b ^ sh transition is given by an 
effective interaction (see e.g., (IH) 



^bsh — 
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SLbRh+{h.c.). (2) 



Using this vertex, Eq. and safely assuming m^ 3> 
we integrate out the massive Higgs boson to obtain the 
effective Lagrangian for the 6 — > s transition with missing 
energy in the final state: 

^b~>sf: = ^CoAimhSLbjiS'^ - C l1 ^,b lv-^ + [h.c). 

(3) 

Leading order Wilson coefficients for the transitions with 
dark matter scalars or neutrinos in the final state are 
given by 



C 



A 



DM 



WwVtlVtb 



9w 



167r2 



xi + 2xt 



+ 



327r2 
3x2 



6xt 



(4) 



{xt-i) 8ixt-iy 



In Xt 



where Xt — mf /M^ . 

We would like to remark at this point that the numer- 
ical value of Cdm is a factor of few larger than C^, 

Cdm 
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if Am^ ~ 0(g^M^ ). This happens despite the fact 
that the effective bsh vertex is suppressed relative to bsZ 
vertex by a small Yukawa coupling ~ m^/vEW- The 
^/vew in lO is compensated by a large coupling of h to 
52, proportional to Xvew, and mh is absorbed into the 
definition of the dimension 6 operator mbSLbRS^ . 

We concentrate on exclusive decay modes with miss- 
ing energy, as these are experimentally more promis- 
ing than inclusive decays and give sensitivity to a large 
range of ms- A limit on the branching ratio has recently 
been reported by BaBar collaboration, Bie+->k+ui> < 
7.0 X lO^'^ at 90% c.l. which improves on a previous 
CLEO hmit [13, but is still far from the SM prediction 
Bi{B Kvv) ~ (3 - 5) X 10-*^ (See, e.g. Q). We use 
the result for C^^sf along with the hadronic form fac- 
tors determined via light-cone sum rule analysis in 
and related to the scalar B ^ K transition in 0|, to 
produce the amplitude of S — > KSS decay. 
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= s ~ {pb — PkY and the form fac- 
tor for B ^ K transition is approximated as /q ~ 
0.3exp{0.63sM^2 - O.mbs^Mg^ -t- O.SQls^M^^}. 

The differential decay width to a, K meson and a pair 
of WIMPs is given by 

xlCl,,h{sf I{s,ms)mliMl 



dV 
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where I(s, ms) reflects the available phase space, 

/(l,ms) = [s2-2J(Af|+Afi)+(M|-A4)2]5[i_4m2/J]i 

From Eq. Q and the prediction for the SM neutrino 
channel, we obtain the total branching ratio for the 
to decay with missing energy in the final state, 



~ 4 X 10"^ + 2.8 X 10-'^K^F{ms) 
Eq. © uses the parametrization of A^m^^, 
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and the available phase space as a function of the un- 
known ms, 



Fims) 



/o(s)^/(s,ms) ds 



/o(s)2/(s,0) ds 



Notice that F(0) = 1 and F{ms) = for ms > ^{ms - 
itt-k) by construction. Similar calculations can be used 
for the decay B K*SS, 

BiB+^K+'+f: ^ 1-3 X 10"^ + 4.3 x IQ- '^n'^F{ms). (10) 

with an analogous form factor. 

For light scalars, ms ~ few 100 MeV, and k 0(1) 
the decay rates with emission of dark matter particles 
are ~ 50 times larger than the decay with neutrinos in 
the final state! This is partly due to a larger amplitude, 
Eq. (O, and partly due to phase space integral that is a 
factor of a few larger for scalars than for neutrinos if ms 
is small. 



ABUNDANCE CALCULATION AND 
COMPARISON WITH EXPERIMENT 

The scalar coupling constant A and the scalar mass 
ms are constrained by the observed abundance of dark 
matter. For low ms, as shown in 4], the acceptable 
value of K is k ^ 0{1). Here we refine the abundance 
calculation for the range < ms < 2.4 GeV in order 
to obtain a more accurate quantitative prediction for k. 
The main parameter that governs the energy density of 
WIMP particles today, which we take to be equal to the 
observed value of Q.DMh? ^ 0.13 0|, is the average of 
their annihilation cross section at the time of freeze-out. 
This cross section multiplied by the relative velocity of 
the annihilating WIMPs is fixed by Q.dm and can be 
conveniently expressed Q as 
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Here F^x denotes the partial rate for the decay, h ^ X, 
for a virtual Higgs, h, with the mass of mf^ = 2Es — "ims- 
Notice that Eq. 1)11(1 contains the same combination 
X^mJ^'^ as (jHJ. The zero-temperature width F^x was ex- 
tensively studied two decades ago in conjunction with 
searches for light Higgs 0, 0, . 

For ms larger than m,r the annihilation to hadrons 
dominates the cross section, which is therefore prone to 
considerable uncertainties. At a given value of ms, we 
can predict F^x within a certain range that reflects these 
uncertainties. With the use of (|ll|l . this prediction trans- 
lates into the upper (A) and the lower (B) bounds on 
K{ms), which we insert into Eq. JH)) and plot the result- 
ing BYB+^K^+f in Figure 13 

In the interval 150 MeV < ms < 350 MeV the an- 
nihilation cross section is dominated by continuum pi- 
ous in the final state, and can be calculated with the 
use of low-energy theorems to good accuracy. Re- 
quiring < 47r allows to determine the lower end of 
the Lee- Weinberg window in our model, TOmin ~ 350 
MeV. In the interval 350 MeV — 650 MeV the strangeness 
threshold opens up, and annihilation into pions via the 
s-channel /o resonance become important. The strength 
of this resonance and its width and position at freeze- 
out temperatures, T ~ (0.05 — 0.1)ms, are uncertain. 
Curve B in this domain of Figure [3 assumes the /o res- 
onance is completely insensitive to thermal effects and 
has the m,inimum width quoted by PDG |2ll | which max- 
imizes Ffcx, whereas curve A corresponds to a complete 
smearing of the /o resonance by thermal effects and much 
lower value of Tf,x- Above ms =1 GeV, curve B takes 
into account the annihilation into hadrons mediated by 
cts{Gfj,u)^ with the two-fold enhancement suggested by 
charmonium decays |l8l | whereas curve A uses the per- 
turbative formula. The charm threshold is treated sim- 
ply by the inclusion of open charm quark production at 
a low threshold (toc ~ 1.2 GeV) in B curve and at a 
high threshold (ms > mjj) in curve A. Both curves in- 
clude the T threshold. There are no tractable ways of 
calculating the cross section in the intermediate region 
650 MeV < ms < 1 GeV. However, there are no particu- 
lar reasons to believe that the annihilation into hadrons 
will be significantly enhanced or suppressed relative to 
the levels in adjacent domains. In this region, we inter- 
polate between high- and low-energy sections of curves 
A and B. Thus, the parameter space consistent with the 
required cosmological abundance of S scalars calculated 
with generous assumptions about strong interaction un- 
certainties is given by the area between the two curves, 
A and B. 

Figure 2 presents the predicted range of Br^+^^+^^i; 
as a function of ms and is the main result of our paper. 
The SM "background" from B Kvv decay is subdom- 
inant everywhere except for the highest kinematically al- 
lowed domain of ms- To compare with experimental re- 
sults [l3. IT3| . we must convert the limit on Br^+^^+^p 
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FIG. 2: Predicted branching ratios for the decay B —> K+ 
missing energy, with current limits from BaBar (I) p^ . 
CLEO (III) O and expected results from BaBar (II). Pa- 
rameter space above curves I and III is excluded. The hori- 
zontal line shows the SM B — > Kvi/ signal. Parameter space 
to the left of the vertical dashed line is also excluded by 
K+ n+f]. 



an efficient probe of dark matter when pair production 
of WIMPs in B meson decays is kinematically allowed. 
In particular, we have shown that the minimal scalar 
model of dark matter with the interaction mediated by 
the Higgs particle predicts observable rates for 
and missing energy. A large portion of the parameter 
space with ms ^ 1 GeV is already excluded by current 
BaBar limits. New experimental data should probe a 
wider range of masses, up to tos ~ 2 GeV. The limits 
obtained in this paper have important implications for 
Higgs searches, as the existence of relatively light scalar 
WIMPs leads to the Higgs decays saturated by invisible 
channel. Given the astrophysical motivations for GeV 
and sub-GeV WIMPs combined with insensitivity of ded- 
icated dark matter searches in this mass range, it is im- 
portant to extend the analysis of 6 ^ s transition with 
missing energy onto other models of light dark matter. 

We thank Misha Voloshin for valuable discussion. This 
research is supported in part by NSERC of Canada and 
PPARC UK. 



to a more appropriate bound on b+^k++$ according 
to the following procedure. We multiply the experimental 
limit of 7.0 x 10~^ by a ratio of two phase space integrals, 
F{ms,Sjnin)/F{ms,Sexp), where Sexp is determined by 
the minimum Kaon momentum considered in the experi- 
mental search, namely 1.5 GeV. This produces an exclu- 
sion curve, nearly parallel to the ms axis at low mg, and 
almost vertical near the experimental kinematic cutoff. 
The current BaBar results (curve I) exclude ms < 430 
MeV, as well as the region 510 MeV< ms < 1.1 GeV, 
and probe the allowed parameter space for dark matter 
up to ms ~ 1.5 GeV. Generalized model with N com- 
ponent dark matter scalar gives A^^-fold increase in the 
branching ratio 0|, and thus greater sensitivity to mg. 

The B factories will soon have larger data samples and 
can extend the search to lower Kaon momenta. The level 
of sensitivity expected from an integrated luminosity C 
of 250 fb~^ and momentum cutoff of 1 GeV is shown 
by curve II, which assumes that the sensitivity scales as 
Cr^l'^ , as suggested by the analysis in 12|. In reality, 
the experimental limit will extend to Kaon momenta be- 
low 1 GeV where the sensitivity will gradually degrade 
due to increasing backgrounds; however, we expect the 
implication of curve II to remain valid, namely that the 
B factories will probe scalar dark matter up to 2 GeV. 

If ms < 150 MeV, the decay K+ n+SS becomes 
possible. The width for this decay can be easily cal- 
culated in a similar fashion to b ^ s transition. The 
concordance of the observed number of events with the 
SM prediction 22] rules out scalars in our model with 
ms < 150 MeV. This exclusion limit is shown by a ver- 
tical line in Figure 3. It is below TOmin of 350 MeV. 

To conclude, we have demonstrated that the 6 — > s 
transitions with missing energy in the final state can be 
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